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LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


A  computer  program  has  been  written  and  has  been  applied  to  LASA  time 
series  data  from  nearly  200  events  in  order  to  obtain  data  for  discrimination  studies. 
In  general  the  previously  estimated  performance  of  M  -m,  ,  spectral  ratio,  and 
complexity  has  been  corroborated.  Some  of  the  specific  results  are  the  following. 
Except  for  one  anomalous  earthquake,  all  shallow  earthquakes  with  m^M.  8  and 
explosions  with  m^^.l  have  been  correctly  identified  using  M^m^.  No  events 
with  m^(5.1  could  be  unequivocally  identified  as  explosions.  However,  the  probab¬ 
ility  that  an  earthquake  would  be  identified  as  such  decreased  to  zero  only  gradually 
as  m^  was  reduced  to  about  4. 0.  A  modification  of  spectral  ratio  has  been  made 
which  introduced  the  option  to  make  no  decision  concerning  the  nature  of  an  event 
when  the  signal -to -noise  ratio  is  not  sufficiently  large.  The  probability  that  no 
decision  will  be  made  is  zero  for  ^>4. 8  and  unity  for  ^<4. 0.  The  probability  of 
correct  identification  is  high  for  events  which  pass  the  signal -to -noise  ratio  test. 

The  period  of  short  period  P-waves  may  be  a  valuable  discriminant  at  low  magnitudes. 
Many  earthquakes  with  magnitudes  below  4. 5  can  be  identified  as  such  using  the 
period  data.  Depth  phase  picks  have  been  made  for  about  60%  of  the  earthquakes 
in  our  population.  About  70%  to  95%  of  these  picks  correspond  to  valid  depth  phases. 
Unfortunately,  depth  phases  were  also  picked  for  three  presumed  explosions. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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I. 


INTRODUCTION 


Short -period  data  from  the  Large  Aperture  Seismic  Array  (LAS A)  has  been 
used  in  the  past  to  develop  and  evaluate  single  station  seismic  discriminants  for  events 

at  teleseismic  distances  from  the  array.  In  1966,  Kelly*  reported  preliminary  LASA 
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results  using  complexity  as  a  discriminant.  Later,  in  1967,  Briscoe  and  Walsh  re¬ 
ported  initial  success  using  a  discriminant  involving  the  frequency  distribution  of  short- 
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period  energy  (spectral  ratio).  More  recently,  Kelly  has  evaluated  these  short -period 

discriminants  using  more  extensive  data. 

The  long-period  data  from  LASA  has  been  used  in  conjunction  with  the  short 

period  data  to  evaluate  the  ability  of  a  single  station  to  identify  events  based  upon  the 

relative  excitation  of  body  phases  and  surface  phases.  (This  is  the  so-called 

Mg  -  m^  criteria. )  Identification  based  on  the  Mg  —  criteria,  but  using  many  tele- 
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seismic  stations,  had  been  previously  reported  by  Liebermann  and  Pomeroy.  ’ 

Earlier  results  relating  to  the  relative  excitation  of  body  and  surface  waves  used  only 
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nearby  stations. 

The  Seismic  Discrimination  Group  at  Lincoln  Laboratory  has  continued  to 
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investigate  the  discriminants  considered  by  Kelly  and  Capon,  et  al.  ’  The  present 
report  is  devoted  primarily  to  further  studies  of  single  station  discriminants  which 
have  been  measured  for  a  large  population  of  almost  200  events. 

It  is  noteworthy  that  the  Lincoln  Laboratory  event  processing  procedures  have 
now  been  considerably  modified  and  standardized.  This  has  greatly  reduced  the  human 
effort  required  to  obtain  standard  discriminants  for  individual  events  and  has  resulted 
in  the  generation  of  a  library  of  events  which  can  conveniently  be  used  to  aid  in  under¬ 
standing  standard  discriminants  or  in  discovering  new  ones.  (The  details  of  the  new 
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event  processing  procedures  can  be  found  in  Section  III. )  All  processing  of  LASA 

waveforms  for  this  report  has  been  accomplished  using  a  single  computer  program 

(SPLP)  designed  to  completely  process  an  event  (both  short-period  and  long-period  data) 

during  a  single  pass  on  a  large  digital  computer. 

The  event  population  discussed  in  this  report  has  been  designed  to  at  least 

partially  rectify  some  of  the  deficiencies  of  the  populations  used  previously.  Previous 

populations,  particularly  for  spectral  ratio,  included  many  large  explosions  and  very 

few  equally  large  earthquakes.  The  Lincoln  Library  of  seismic  data  has  been  searched 

for  large  earthquakes  which  can  be  compared  with  the  large  explosions.  The  population 
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®f  presumed  explosions  has  been  increased  nearly  100%  above  that  used  by  Kelly. 

Whenever  possible,  both  long-period  and  short -period  waveforms  are  now 

processed  for  each  event.  Thus  this  data  can  be  considered  jointly  without  great  dif- 
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ficulty.  The  long-period  and  short-period  data  used  in  earlier  reports  ’  ’  was  ob¬ 
tained  somewhat  independently  making  it  difficult  to  obtain  a  common  data  base.  In 
addition,  the  procedures  for  data  reduction  which  are  now  in  use  are  uniformly  applied 
to  all  events.  Thus  nuisance  effects  due  to  minor  variations  in  the  treatment  of  indi¬ 
vidual  events  have  been  minimized  if  not  eliminated. 

All  measured  diagnostics  as  well  as  short -period  LASA  beams  and  other  data 
for  all  events  in  our  population  are  stored  on  a  single  composite  digital  magnetic  tape. 
New  events  can  be  added  as  desired.  This  digital  tape  storage  of  data  has  greatly  aided 
in  the  study  of  standard  discriminants.  In  addition,  the  data  on  the  composite  tape  is 
currently  being  used  to  investigate  possible  multivariate  use  of  the  data  from  a  single 
site. 
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The  same  program  (SPLP)  which  has  been  used  to  generate  the  composite  tape 
of  LASA  data  can  be  applied  to  the  short-period  data  currently  being  recorded  at 
NORSAR.  Thus  a  composite  tape  of  data  recorded  in  Norway  can  be  generated.  This 
is  currently  in  progress  for  that  subset  of  events  in  our  data  base  which  might  have 
been  sensed  in  Norway  and  for  which  recordings  were  made.  Ultimately  it  is  hoped 
that  the  LASA  and  NORSAR  composite  tapes  can  be  used  for  a  joint  discrimination  ex¬ 
periment  in  which  data  processing  at  both  sites  is  as  uniform  as  is  feasible. 
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II. 


DATA  BASE 


The  earthquake  population  considered  in  this  report  consists  primarily  of  events 
believed  to  be  in  or  very  close  to  the  Sino-Soviet  bloc  countries.  In  addition,  the  popu¬ 
lation  includes  four  earthquakes  in  North  Africa  (Algeria,  Tunisia),  10  earthquakes  in 
the  Mediterranean  and  Aegean  Sea  area,  six  earthquakes  in  the  area  of  Iran,  Iraq  and 
western  Thrkey,  and  six  earthquakes  from  the  Laptev  Sea  and  Arctic  Ocean.  Altogether 
the  earthquake  population  consists  of  156  events.  Of  these  about  40%  are  in  the  Kurile  - 
Kamchatka  area. 

The  population  of  35  presumed  underground  nuclear  explosions  includes  events 
from  Sahara  and  Amchitka  (Longshot)  in  addition  to  events  from  Soviet  test  sites. 

Our  earthquake  population  is  not  representative  of  the  natural  seismicity  of  the 
earth.  For  example,  LASA  station  bulletins  were  initially  searched  starting  1  January 
1968  for  events  in  or  near  the  Sino-Soviet  bloc  countries.  This  population  began  to 
be  biased  towards  events  with  magnitudes  from  4.  0  to  5.  0  on  the  Richter  scale  and 
towards  events  from  the  Kurile -Kamchatka  area.  Events  were  then  added  only  if  they 
were  large  or  from  outside  the  Kurile -Kamchatka  area.  It  was  felt  that  the  population 
was  still  deficient  in  large  magnitude  events  so  the  entire  Lincoln  Laboratory  tape  li¬ 
brary  was  searched  for  large  events  in  the  Sino-Soviet  area  which  could  be  added. 

Some  of  these  large  earthquakes,  as  well  as  many  of  the  nuclear  explosions, were  part 
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of  the  data  base  for  the  experiments  previously  reported  by  Kelly  and  by  Capon,  et  al. 
Figure  1  shows  a  cumulative  histogram  of  the  magnitude  of  the  earthquake  population. 

A  line  with  the  same  slope  as  the  natural  seismicity  curve  has  been  drawn  on  the  figure 
to  emphasize  that  our  population  is  not  representative  of  natural  seismicity. 


Only  about  55%  of  the  earthquakes  and  about  70%  of  the  presumed  underground 
explosions  included  in  our  data  have  been  reported  by  the  U.  S.  Coast  and  Geodetic 
Survey  (USCGS).  This  does  not  imply  anything  about  the  relative  detection  levels  of 
LASA  and  USCGS.  It  does  define  the  size  of  the  population  which  we  have  used  most 
extensively  for  checking  the  ability  of  a  single  large  array  to  determine  depth. 

Our  population  of  events  includes  many  which  are  anomalous  in  one  or  more 
aspects.  For  example,  consider  Fig.  2  which  is  a  plot  of  the  LASA  body-wave  magni¬ 
tudes  against  the  USCGS  magnitudes  for  both  earthquakes  and  explosions.  Events  more 
than  99°  from  LASA  or  with  USCGS  depths  greater  than  100  km  have  not  been  shown. 
There  is  general  agreement  between  the  two  magnitudes  with  LASA  tending  to  be 
slightly  higher.  The  scatter  of  the  main  body  of  data  is  generally  less  than  ±  0.  5  mag¬ 
nitude  units.  A  few  points  do  scatter  to  magnitude  discrepancies  as  large  as  1. 0  mag¬ 
nitude  units.  This  data  is  included  here  to  underline  the  fact  that  even  using  a  large 
array  to  measure  an  event  property  as  basic  as  magnitude  will  result  in  some  scatter 
and  occasionally  very  anomalous  measurements.  No  attempt  has  been  made  to  exclude 
events  from  our  data  base  because  of  some  anomalous  behavior  or  properties.  Indi¬ 
vidual  events  are  discussed  in  other  sections  only  when  their  anomalous  behavior 
relates  to  discrimination.  Thus  an  event  with  a  large  LASA-USCGS  magnitude  discrep¬ 
ancy  would  not  be  specifically  discussed  unless  the  discrepancy  resulted  in  discrimin¬ 
ation  difficulties. 
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III. 


ROUTINE  EVENT  PROCESSING 


Short-period  and/or  long-period  seismograms  for  all  events  in  our  data  base 
are  recorded  on  standard  format  LASA  digital  tapes.  As  mentioned  in  Section  I,  a 
single  computer  program  (Short-Period  Long-Period  or  SPLP)  has  been  written  which 
normally  allows  us  to  complete  all  routine  analyses  of  LASA  data  for  an  event  in  a  sin¬ 
gle  computer  run  with  the  LASA  digital  tapes  for  the  event.  In  this  Section  we  describe 
the  operational  characteristics  of  this  program  and  the  interaction  of  an  analyst  with 
the  program. 

SPLP  will  accept  either  Lo-Rate  or  Hi -Rate  tape  from  LASA  as  input.  If  Lo- 
Rate  tape  is  supplied,  then  both  long-period  and  short -period  analysis  is  done.  The 
short -period  data  in  this  case  are  the  21  straight  sum  traces  recorded  on  the  tape.  If 
Hi-Rate  tape  is  supplied,  the  program  operates  with  the  twenty-one  500-ft  seismom¬ 
eters  at  the  center  of  each  subarray  and  processes  only  short -period  data.  The  pro¬ 
gram  requires  about  20  minutes  of  IBM  360/67  time  to  completely  process  an  event. 
This  is  reduced  to  10  —  15  minutes  if  only  short-period  data  are  available.  It  is  likely 
that  by  restructuring  and  rewriting  the  program  the  processing  time  could  be  reduced. 
An  absolute  limitation  to  time  reduction  is  imposed  by  the  approximately  eight  minutes 
required  to  read  through  a  complete  input  tape.  SPLP  also  accepts  the  short-period 
data  generated  by  the  limited  array  in  the  area  of  Oyer,  Norway. 

Figure  3  is  a  flow  chart  indicating  all  of  the  major  operations  done  by  SPLP. 
Except  for  the  entry  of  seismograms  into  core,  no  input  or  output  has  been  shown.  All 
measured  or  computed  event  parameters  as  well  as  beam  waveforms  and  spectral 
curves  are  output  graphically  using  an  SC  4020  hard  copy  display  and  are  also  stored 
on  a  composite  digital  magnetic  tape.  Although  many  options  are  available,  only  a 
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small  number  of  input  parameters  is  changed  for  the  routine  processing  of  events.  In 
particular,  the  program  requires  as  an  input  rough  estimates  of  P  arrival  time,  azi¬ 
muth,  and  horizontal  phase  velocity  for  the  event.  Arrival  times  within  20  seconds 
are  quite  sufficient,  and  the  estimated  azimuth  (3)  and  horizontal  phase  velocity  (v) 
supplied  to  the  program  can  be  quite  crude.  The  arrival  times,  Q,  and  v  obtained 
from  LASA  station  bulletins  have  proved  to  be  more  than  sufficiently  accurate.  If  the 
quality  of  0  and  v  are  thought  to  be  particularly  poor,  a  larger  grid  of  SP  beams  can 
be  requested. 

The  program  can  be  requested  to  use  any  combination  of  long-period  and 
short -period  instruments  to  perform  an  analysis.  Thus,  for  example,  when  it  is 
known  that  a  particular  instrument  is  not  operational,  the  program  is  requested  to 
omit  it  from  consideration.  However,  the  program  was  designed  to  operate  without 
the  need  to  visually  check  every  waveform  to  be  used.  Some  simple  checks  are  in¬ 
cluded  to  detect  and  reject  bad  data.  The  peak  value  and  average  absolute  value  for 
each  channel  of  data  are  computed  over  the  time  interval  of  interest.  Any  channel  for 
which  either  the  peak  or  the  average  differs  by  more  than  a  factor  of  five  from  the 
median  channel  is  removed  from  consideration.  This  procedure  has  been  quite  suc¬ 
cessful  in  removing  anomalous  channels.  More  sophisticated  tests  of  data,  with  pos¬ 
sible  data  correction  for  isolated  glitches  in  the  waveform,  might  have  saved  a  few  of 
the  channels  rejected  by  the  program  and  eliminated  the  cases  in  which  intervention  by 
an  analyst  was  required  to  select  satisfactory  channels. 

Bandpass  filtering  of  all  long-period  and  short -period  data  has  been 
accomplished  in  the  frequency  domain  as  described  in  Appendix  A.  The  0. 6  —  2. 0  Hz 
and  0. 025  —  0. 055  Hz  passbands  indicated  on  Fig.  3  have  been  used  for  all  events. 
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These  passbands  evolved  from  earlier  experience  with  LASA  data  and  test  runs  with 
SPLP.  The  bands  can  be  changed  by  data  cards  submitted  with  SPLP.  This  is  convenient 
since  for  example  initial  experience  with  Norway  data  indicates  that  a  0.  6  —  3. 5  Hz 
passband  yields  much  more  satisfactory  operation  than  does  0. 6  —  2. 0  Hz,  which  often 
seriously  degrades  the  signal.  One  might  adjust  the  passband  for  each  event  but  that 
would  require  significant  pre-run  analysis  of  the  event. 

Let  vQN  and  VqE  be  the  north  and  east  components  of  velocity  corresponding  to 
the  3  and  v  input  to  SPLP  for  an  event.  The  program  forms  a  filtered  short -period  beam 
for  an  event  with  horizontal  phase  velocity  components,  Vq^  and  Vqe*  LASA  station 
corrections  are  added  to  the  plane  wave  moveouts  when  forming  the  beam.  A  number 
of  similar  beams  are  formed  about  this  central  beam.  Beams  are  formed,  including 
station  corrections,  with  velocity  components 


and 


VN  =  £o 


^ON  + 


(kQN  +  mAk)2  +  (kQE  +  nAk)2 


v,-,  =  f 
E  o 


k0E+  nAk 


(kQN  +  mAk)2  +  (kQE  +  nAk)2 


where  f  is  a  nominal  event  frequency  read  into  the  program,  Ak  is  a  wavenumber 
increment  read  into  the  program,  n  and  m  are  positive  or  negative  integers  with  limits 
read  into  the  program. 
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The  grid  of  beams  formed  is  uniform  in  wavenumber.  It  is  not  uniform  in 
location  or  velocity.  The  density  of  beams  in  wavenumber  is  fixed  by  Ak  which  is  an 
input  parameter.  Most  runs  with  LASA  data  have  been  made  with  Ak  =  0.  00125  c/km. 
The  LASA  3  db  beam  width  is  about  0. 005  c/km.  Thus  several  beams  will  be  formed 


within  the  main  LASA  beam.  If  an  event  is  85°  distant  from  LASA,  the  choice 

Ak  =  0. 00125  results  in  beams  separated  by  roughly  2°  on  the  earth's  surface.  This  is 


typical  of  the  separation  between  beams. 

Event  location  is  done  by  choosing  that  filtered  beam  with  the  largest  peak  to 
peak  excursion.  Zero  depth  is  assumed  and  the  velocity  for  that  beam  is  converted  to 
bearing  and  distance  from  LASA.  The  resulting  epicenter  is  not  the  best  obtainable 
using  LASA.  The  grid  density  and  location  procedure  have  been  chosen  to  yield  a  lo¬ 
cation  sufficient  for  our  discrimination  studies  not  to  yield  a  best  epicenter.  Our  cur¬ 
rent  objective  has  not  been  to  investigate  discrimination  based  upon  accurate  epicenter 
determination. 

The  program  uses  the  best  filtered  beam  to  estimate  body-wave  magnitude.  The 
largest  peak-to-peak  excursion  in  the  waveform  and  its  period  are  measured.  From 
this  and  distance  the  LASA  magnitude  is  computed.  ^ 

Complexity  is  measured  by  SPLP  using  the  best  filtered  beam.  The  definition 
of  complexity  used  is 


35 
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C  =  j  |x(t)|  dt  /  J*  |x(t)|  dt 
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where  x(t)  is  the  beam,  the  origin  of  time  is  at  the  event  onset,  and  in  practice  the 
integrals  are  actually  sums  over  data  samples.  The  program  has  no  algorithm  for 
determining  onset  time.  Only  the  time  of  the  peak  excursion  of  the  beam  is  known. 

Thus  the  program  computes  several  complexity  values  assuming  various  times  from 
onset  to  peak  value.  An  analyst  has  ultimate  responsibility,  guided  by  SC  4020  plots 
with  time  intervals  indicated,  for  choosing  which  is  the  correct  measurement.  In 
some  cases,  for  example  a  very  emergent  event,  no  complexity  values  are  valid.  If 
this  is  the  case,  the  correct  value  is  computed  during  a  final  many  event  editing  run  of 
the  SPLP  composite  tape.  The  output  tape  from  this  editing  run  is  described  at  the  end 
of  this  Section  and  in  Appendix  B. 

A  very  preliminary  check  for  possible  depth  phases  is  done  by  SPLP.  A  short 
section  of  beam,  typically  five  seconds,  around  the  peak  time  is  zeroed  out.  The  re¬ 
maining  waveform  is  scanned  to  find  the  location  of  the  largest  peak  which  is  greater 
than  0.  25  of  the  peak  in  the  blanked  out  interval.  Such  a  peak  is  indicated  on  the  SC 
4020  output.  Whenever  possible,  USCGS  data  is  obtained  for  each  event.  If  this  data 
indicates  a  possible  depth  phase  outside  of  the  102. 8  second  short -period  window  used 
by  SPLP  a  second  run  is  undertaken  explicitly  to  detect  any  possible  depth  pulse.  Ulti¬ 
mate  responsibility  for  all  phase  picking  is  given  to  the  analyst.  His  decisions  about 
possible  depth  phases,  if  he  is  at  variance  with  the  program  pick,  are  included  during 
the  final  editing  run. 

The  epicenter  determined  from  filtered  data  is  used  to  form  an  unfiltered  short 
period  beam  waveform.  This  beam  waveform  is  used  to  obtain  spectra  and  to  measure 
spectral  ratio.  Spectra  have  been  computed  for  10  and  20  second  intervals  including 
the  event.  The  spectra  are  computed  by  discrete  Fourier  transforming  the  entire  beam. 
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normally  102. 8  seconds  long,  with  all  data  outside  the  10  or  20  second  interval  set  to 
zero.  The  data  is  tapered  linearly  to  zero  over  1.0  seconds  at  the  ends  of  the  intervals. 
Just  as  is  the  case  for  complexity  the  onset  time  is  not  known.  Thus  several  shifted 
intervals  were  used  and  the  analyst  picked  the  appropriate  one.  If  none  were  satisfac¬ 
tory,  this  was  corrected  during  the  final  composite  tape  editing  run. 

Spectral  ratio  is  computed  for  both  the  10  second  and  20  second  intervals.  The 

3  4 

definition  has  been  changed  from  that  used  previously  by  Briscoe  and  Kelly.  Specif¬ 
ically,  SPLP  computes 

1.95  0.85 

R  =  J  |  X  (f )  |  df  /  J‘  |  X  (f )  |  df 

1.45  0.35 

where  |X(f)|  is  the  magnitude  of  the  10  or  20  second  interval  spectra  and  in  practice 
the  integrals  are  replaced  with  sums  over  frequency  components. 

Travel  time  tables  for  Rayleigh  as  well  as  P  waves  are  included  in  SPLP.  Thus 
once  all  the  short -period  data  is  processed,  it  automatically  moves  to  the  appropriate 
time  interval  and  begins  processing  long-period  data. 

Beamforming  of  long-period  data  is  done  with  no  station  corrections  and  no 
correction  for  the  slightly  different  velocities  of  the  different  frequency  components 
of  the  long-period  events.  Both  of  these  refinements  would  result  in  at  best  a  trivial 
improvement  of  the  beams.  The  nominal  Rayleigh  wave  velocity  used  for  steering  has 
been  determined  by  measuring  events.  The  surface  waves  are  assumed  to  arrive  from 
the  same  azimuth  as  the  P  waves. 

Unlike  earlier  experiments^  all  chirp  filtering  is  done  in  the  frequency  domain 
as  is  described  in  Appendix  A.  Chirp  filters  with  impulse  response  durations  ranging 
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from  200  seconds  to  800  seconds  in  increments  of  50  seconds  are  applied  to  the  vertical 
long-period  beam.  Matched  filtering  such  as  that  used  by  Alexander^  has  not  been 
used.  The  chirp  filtered  beam  with  the  largest  zero  to  peak  excursion  is  saved  to 
compute  a  number  which  should  differ  from  surface  wave  magnitude  by  only  an  additive 
constant. 

Specifically,  SPLP  computes 

L  =  log1Q  A  + 1.66  log1Q  A 

where  A  is  the  amplitude  of  the  largest  chirp  waveform  and  A  is  epicentral  distance. 

12 

The  relation  of  L  to  M  ,  the  surface  wave  magnitude,  is  discussed  in  Section  IV. 

s 

The  analyst,  utilizing  plots  of  all  long-period  beams  including  chirp  filtered 
radial  and  vertical  beams,  has  the  ultimate  judgments  to  make  about  the  long-period 
data.  He  must  decide  if  (1)  the  program  has  picked  the  event  properly,  (2)  no  signal 
is  actually  visible  above  the  background  noise,  or  (3)  some  other  event  has  arrived  at 
about  the  same  time  and  obscured  the  event  of  interest. 

Much  of  the  data  recorded  on  the  SPLP  composite  tape  is  redundant,  unwanted, 
obsoleted  by  later  data  on  the  same  tape,  or  should  be  modified  or  corrected  by  inputs 
from  an  analyst.  This  is  done  simultaneously  for  all  events  on  the  composite  tape 
using  edit  programs  designed  for  the  purpose.  The  output  of  the  edit  program  is  an 
800  bpi  9-track  tape  with  a  single  logical  record  for  each  successfully  processed  event 
on  the  SPLP  composite  tape.  Each  logical  record  is  written  by  an  IBM  360  Fortran  IV 
write  statement.  The  record  format  is  specified  as  variable  with  a  blocksize  of  2688 
bytes  and  a  record  length  of  2684  bytes.  It  is  this  tape  which  we  have  used  for  our  dis¬ 
crimination  studies.  Appendix  B  shows  the  structure  of  each  logical  record.  Copies  of 
this  edited  tape  are  available  to  qualified  requesters  of  LASA  data. 
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IV. 


CALIBRATION  FOR  SURFACE  WAVE  MAGNITUDE 


The  SPLP  program  does  not  directly  generate  a  surface  wave  magnitude.  It 
does  generate  a  surface  wave  factor,  L,  which  was  defined  in  Section  III.  For  A 
between  15°  and  130°  Gutenberg^  gives  surface  wave  magnitude  as 

Mg  =  log^Q  A'  +  1.  66  log-^Q  A  +  1.  82 

where  A'  is  the  horizontal  component  of  the  maximum  ground  movement  in  microns 
during  surface  waves  having  periods  of  about  20  seconds.  This  expression  for  M 

s 

ignores  station  effects  and  the  depth  of  events.  The  surface  factor  and  surface  wave 
magnitude  are  approximately  related  by 

M  =  L+ 1.7 
s 

This  relation  was  determined  by  the  calibration  procedures  described  below  and  has 

been  used  throughout  this  report  to  obtain  M  values  for  events. 

s 

One  approach  is  to  assume  that  the  horizontal  and  vertical  components  of  the  LP 
Rayleigh  wavetrain  will  have  approximately  the  same  amplitude.  Assume  also  that  the 
peak  of  the  best  chirp  filtered  wavetrain  is  B  times  the  peak  of  the  input  wavetrain. 

That  is,  we  assume  that  BA’  =  A,  where  A  is  the  amplitude  of  the  chirp  output.  It  then 
follows  that 


Mg  =  L  +  1.  82  —  logjQ  B 

We  have  measured  B  for  a  collection  of  33  earthquakes  and  10  bombs  which  had 
clear  surface  waves  before  chirp  filtering.  The  value  of  B  obtained  for  earthquakes  was 
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1. 48  and  for  the  explosions,  1. 68.  The  overall  average  value  of  B  was  thus  1.  53.  Using 
this  value  of  B  gives  Mg  =  L  +  1.  63. 

Another  calibration  approach  has  been  to  compare  L  with  M  values  independently 

s 

reported  for  the  same  events.  This  has  been  done  for  10  earthquakes  and  seven  explo¬ 
sions  in  our  data  base  which  have  had  surface  wave  magnitudes  independently  measured. 

The  events  for  comparison  had  M  determined  previously  using  either  unprocessed 

s 

5  6 

seismograms  or  time -domain  chirp  filtering  methods.  '  This  comparison  gave 

M  =  L  +  1.72. 
s 

The  relationship  Mg  =  L  +  1.  7  has  been  accepted  as  a  reasonable  compromise 

between  those  obtained  by  the  above  calibration  procedures. 

Figure  4  is  a  plot  of  M  =  L  +  1.  7  versus  LASA  body -wave  magnitude  for 

s 

earthquakes  in  our  population.  The  plot  includes  a  number  of  events  never  reported  by 
USCGS,  but  all  events  located  deeper  than  100  km  by  USCGS  have  been  excluded.  The 
empirical  relationship  Mg  =  1. 59  m^  —  3.  97,  determined  by  Gutenberg  and  Richter^ 
has  been  drawn  on  the  figure.  This  line  is  a  reasonable  average  line  for  the  data  shown. 
This  plot,  which  is  in  fact  a  somewhat  different  method  of  calibration,  verifies  the 
validity  of  the  relationship  of  Mg  and  L. 
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V. 


SURFACE  WAVE  DETECTION  PROBABILITY 


Capon,  et  al,'^  have  previously  considered  the  same  geographical  area  as  we  are 
considering  and  reported  that  surface  waves  were  detected  by  LASA  for  all  earthquakes 
with  body-wave  magnitudes  greater  than  about  4.  9.  Our  data  has  substantiated  those 
results  by  achieving  a  100%  detection  threshold  of  4.  8.  In  addition,  we  have  considered 
our  data  in  sufficient  detail  to  obtain  estimates  of  the  probability  of  detection  of  surface 
waves  for  each  value  of  body -wave  magnitude. 

Figure  5  shows  three  incremental  histograms  for  three  sets  of  events  in  our 
population.  The  top  is  for  all  earthquakes  for  which  we  detected  surface  waves.  The 
bottom  plot  is  for  all  earthquakes  which  had  surface  waves  obscured  by  the  simultaneous 
arrival  of  surface  wave  trains  from  comparable  or  larger  events.  The  center  plot  is 
for  earthquakes  with  surface  waves  obscured  by  the  continuous  background  seismic 
noise.  The  magnitude  boxes  used  to  construct  these  histograms  were  each  0. 1  body- 
wave  magnitude  unit  wide.  No  histogram  is  shown  of  the  events  for  which  equipment  or 
tape  problems  precluded  the  detection  of  surface  waves. 

Figure  6  was  constructed  from  the  upper  two  histograms  on  Fig.  5.  The  dots 
represent  the  incremental  detection  probability  computed  by  taking  the  ratio  of  detected 
events  to  total  events  in  each  magnitude  range.  Events  deeper  than  100  km  were  exclu¬ 
ded  in  constructing  Fig.  6.  The  scatter  of  the  data  is  large.  The  100%  detection  at 
m^  =  4. 8  is  probably  optimistic  although  it  appears  reasonable  to  accept  a  value  as 
large  as  90%. 

The  detection  probability  can  be  estimated  in  a  different  way  using  measured 
noise  properties.  The  approach  is  much  more  indirect.  Nevertheless,  the  results  are 
not  inconsistent  with  the  direct  estimates. 
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Indirect  estimates  were  obtained  as  follows.  First,  all  events  with  no  detected 

surface  waves  were  used  to  construct  Fig.  7.  The  figure  was  constructed  by  using  the 

largest  noise  burst  and  the  distance  to  die  event  to  compute  M  values.  Although  the 

s 

resulting  distribution  of  M  values  may  be  biased  high,  Fig.  7  is  interpreted  to  be  die 

s 

distribution  of  detectable  surface  wave  magnitudes  if  detection  is  always  possible  at 

0.  0  db  signal -to-noise  ratio.  Thus  for  example,  if  detection  is  always  possible  at 

0.  0  db,  then  approximately  50%  of  all  events  with  surface  wave  magnitude  2.  9  will  be 

detected  since  the  noise  will  be  less  than  that  equivalent  M  value  50%  of  the  time.  If 

s 

a  different  SNR  is  required  for  detection,  then  the  curve  of  Fig.  7  must  be  shifted  to 

obtain  the  probability  of  detection  for  a  given  M  value.  For  example,  if  3  db  of  SNR 

s 

is  required,  the  curve  must  be  moved  to  the  right  by  3/20  surface-wave  magnitude 
units.  In  this  case  the  detection  probability  for  events  with  Mg  =  2.  9  is  reduced  from 
0.  5  to  about  0.  25.  Once  the  SNR  required  for  detection  is  established  it  is  now  easy  to 
obtain  the  probability  of  detection  of  surface  waves  as  a  function  of  m,  rather  than  M^. 
Figure  8  shows  such  curves  for  3.  0  db  and  6.  0  db  SNR  required  for  detection.  Those 
curves  are  simply  generated  from  the  appropriately  shifted  curves  of  Fig.  7  by  assuming 
the  Gutenberg -Richter  relationship  Mg  =  1.  59  m^  —  3.  97.  ^ 

There  is  fair  agreement  between  the  detection  probabilities  shown  in  Fig.  8 
(3  db  curve)  and  in  Fig.  6.  Both  curves  indicate  that  the  incremental  detection  proba¬ 
bility  for  earthquakes  is  0.  5  for  about  equal  to  4.  5  and  increases  to  nearly  1.  0  at 
m^  =  4.  8.  The  3.  0  db  SNR  required  for  detection  is  not  unreasonable  for  an  analyst 
who  can  make  use  of  time  windows  and  his  ability  to  recognize  waveform  characteristics. 

The  preceding  discussion  of  detection  probabilities  did  not  consider  events 
obscured  by  interfering  events.  It  is  clear  from  Fig.  5  that  interference  can  be  a 
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nontrivial  problem.  However,  since  our  data  is  not  particularly  representative  of 
natural  seismicity,  we  have  not  attempted  to  consider  the  problem  in  great  detail. 
Our  data  simply  sets  the  order  of  magnitude  of  the  problem. 
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VI.  M  -  m,  DISCRIMINATION 
s  b 

The  relative  excitation  of  surface  waves  and  body  waves  by  earthquakes  and 
underground  explosions  has  previously  been  established  as  a  powerful  discriminant.  The 
data  presented  in  this  Section  essentially  corroborates  the  usefulness  of  the  M  -  m, 

S  D 

discriminant. 

Figure  9  shows  a  scatter  diagram  of  Mg  vs  m^  for  69  earthquakes  with  detected 
surface  waves  and  10  presumed  explosions  with  detected  surface  waves.  Both  and  m^ 
are  the  values  measured  using  SPLP.  There  is  a  good  separation  of  the  two  populations 
except  for  a  few  earthquakes  which  appear  to  have  low  surface -wave  magnitudes. 

Figure  10  shows  the  scatter  diagram  after  removing  the  five  events  in  the  population  of 
Fig.  9  which  were  reported  at  depths  greater  than  100  km  by  USCGS.  In  addition,  M 

o 

upper  bounds  for  five  presumed  explosions  with  surface  waves  obscured  by  noise  have 

been  shown  on  the  plot.  The  separation  between  bombs  and  earthquakes  is  now  clearly 

better.  A  comparison  of  the  two  figures  shows  that  all  events  removed  by  the  depth 

filter  had  relatively  small  surface  waves  resulting  in  small  values  of  M  .  This  is  to  be 

s 

expected  since  deep  events  should  be  less  efficient  producers  of  surface  waves  than 
events  at  normal  depths. 

One  earthquake  remains  entirely  within  the  explosion  region  of  Fig.  10.  That 
event  is  located  just  west  of  the  Kurile  Islands.  The  P-wave  arrival  at  LASA  was  at 
16:14:50  on  23  December  1967.  It  has  been  reported  and  located  by  USCGS  as  well  as 
by  LASA.  The  event  is  in  fact  not  a  real  discrimination  problem  since  it  is  located  well 
off  of  any  land  by  USCGS  and  there  is  a  clear  arrival  at  LASA  which,  if  it  is  assumed  to 
be  pP,  locates  the  event  at  about  33  km  depth.  The  depth  is  corroborated  by  the  26  km 
unrestrained  depth  obtained  by  USCGS.  In  addition,  the  first  motion  recorded  at  LASA 
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is  clearly  dilatational  (Event  E  on  Fig.  26.)  Despite  all  of  this  the  event  is  a 

problem  since  such  a  shallow  event  should  more  effectively  generate  surface  waves. 

The  USCGS  has  assigned  a  body-wave  magnitude  of  5. 1  to  the  problem  event 

whereas  LASA  has  assigned  it  5.  7.  Since  the  USCGS  m^  would  make  the  event  more 

like  an  earthquake  it  was  decided  to  plot  the  LASA  Mg  against  the  USCGS  m^  for  all 

events  on  Fig.  10  which  were  also  reported  by  USCGS.  Figure  11  is  that  plot.  Although 

the  problem  event  still  has  an  unusually  low  value  of  M  the  separation  between  earth- 

s 

quakes  and  explosions  is  now  complete.  It  should  probably  be  noted  that  the  m^  =  5.0 

explosion  had  marginally  detected  Rayleigh  waves.  The  measured  M  for  that  event 

s 

may  be  slightly  high  due  to  noise. 

A  comparison  of  Figs.  10  and  11  suggests  that  the  m^  values  obtained  from  a 

worldwide  system  (USCGS  in  this  case)  is  more  stable  and  effective  for  discrimination 

than  that  obtained  from  a  single  large  array  station.  Various  arguments  revolving 

around  P-wave  radiation  patterns  can  be  employed  to  substantiate  such  claims.  The 

experimental  data  at  hand  is  certainly  not  conclusive.  For  example,  since  USCGS  did 

not  report  all  events  in  our  population.  Fig.  11  contains  many  fewer  points  than  Fig.  10. 

The  value  of  M  measured  at  a  site  is  also  a  function  of  the  radiation  pattern  of 
s 

the  event.  This  might  also  be  affecting  the  discrimination  ability  of  LASA  based  upon 

M  —  m.  . 
s  b 

5  6 

It  is  generally  believed0’  that  the  Mg  —  m^  criteria  should  be  regionalized. 

That  is,  only  explosions  and  earthquakes  from  the  same  tectonic  region  should  be  com¬ 
pared.  This  is  not  contradicted  by  our  data  but  neither  is  it  very  strongly  substantiated. 

For  example,  the  percentages  of  normal  depth  events  above  and  below  the  line  M  = 

s 

1.  59  —  3. 97  have  been  computed.  Using  the  LASA  gives  57%  below  the  line  for 
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Kurile -Kamchatka  earthquakes  and  only  33%  for  the  remaining  earthquakes  in  our 
population.  If  the  USCGS  m^  is  used,  the  figures  become  44%  and  30%,  respectively. 
This  tends  to  substantiate  the  assertion  that  events  from  the  Kurile -Kamchatka  area 
often  have  small  surface  waves.  The  strongest  support  for  regionalization  in  our 
data  may  be  the  23  December  1967  problem  event  which  is  from  the  Kurile  Islands 
region.  None  of  the  presumed  explosions  with  detected  surface  waves  are  from  that 
region. 

It  appears  that  the  Mg-m^  criterion  can  be  successfully  applied  to  most 
earthquakes  with  mj_>  4. 8  since  Mg  can  be  measured  for  such  events.  As  m^  is 
reduced  to  about  4.  0,  the  probability  of  identification  is  reduced  to  zero  as  the  prob¬ 
ability  of  surface  wave  detection  goes  to  zero.  All  presumed  explosions  in  our  data 
base  with  m^>  5. 1 ,  exclusive  of  those  subjected  to  interfering  events  or  with  no  tape 
available,  have  been  correctly  identified  using  M  -m,  .  It  has  not  been  possible  to 

apply  M  -m,  to  any  of  the  smaller  presumed  explosions.  It  would  thus  appear  that 
S  D 

many  earthquakes  can  be  identified  as  such  using  M  -m,  at  lower  body  wave  magni - 

S  D 

tudes  than  those  at  which  explosions  can  be  identified  using  M  -m. . 
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VII.  BASIC  SPECTRAL  RATIO  DISCRIMINATION 


As  mentioned  in  Section  III  spectral  ratio  has  been  calculated  using  a  10-second 
interval  and  a  20-second  interval  of  data  for  each  event.  The  two  spectral  ratios 
were  approximately  equal  in  most  cases.  In  general,  we  have  found  no  significantly 
better  discrimination  using  one  rather  than  the  other  of  these  two  spectral  ratios. 

Since  the  10-second  measure  will  tend  to  have  a  slight  signal -to-noise  advantage  and 
will  less  often  include  clear  depth  phases  in  the  interval  we  have  tended  to  use  it 
rather  than  the  20-second  measure.  In  the  remainder  of  this  section  only  the  10-second 
spectral  ratio  is  considered.  A  comparison  of  10  and  20  second  results  is  given  in 
Section  VIII  after  events  have  been  removed  from  the  population  by  virtue  of  unsatis¬ 
factory  signal-to-noise  ratios. 

Figure  12  is  a  scatter  diagram  of  spectral  ratio  vs  the  LASA  for  all  under¬ 
ground  explosions  and  all  earthquakes  with  distances  less  than  99°.  The  decision  line 
shown  on  the  figure  is  an  exponential  and  is  the  geometrical  mean  of  the  two  indicated 
population  trends.  None  of  these  lines  has  been  formally  optimized  but  were  obtained 
as  described  in  the  following  paragraph. 

The  spectral  ratio  of  an  event  is  the  ratio  of  two  quantities  derived  from  the 
spectrum  of  the  event.  Let  V^j  and  V^  be  these  quantities  obtained  from  the  1.  45  — 

1.  95  Hz  (high)  band  and  the  0. 35  —  0.  85  Hz  (low)  band,  respectively.  Figures  13  and 
14  show  V^j  and  V^  scattered  against  LASA  m^  on  semilog  paper.  These  figures  were 
used  to  obtain  the  decision  lines  and  trend  lines  shown  on  Fig.  12.  The  trends  of  the 
data  on  Figs.  13  and  14  appear  to  be  well  represented  by  straight  lines.  The  straight 
lines  drawn  on  the  figures  were  selected  visually  to  obtain  a  good  fit  to  the  data.  A 
nominal  spectral  ratio  as  a  function  of  magnitude  was  obtained  for  each  type  of  event  by 


21 


using  the  appropriate  trend  lines  on  Figs.  13  and  14.  These  nominal  spectral  ratios 
are  the  trend  lines  shown  on  Fig.  12.  Note  that  spectral  ratio  discrimination  is  equiv¬ 
alent  to  using  log  V^j  —  log  as  a  discriminant.  The  decision  line  shown  on  Fig.  12 
is  equivalent  to  using  the  arithmetic  means  of  the  log  of  spectral  ratio  trends  as  a 
decision  line  and  the  log  of  spectral  ratio  as  a  discriminant. 

Figure  15  shows  a  plot  of  spectral  ratio  vs  LASA  m^  for  the  same  events  as 
Fig.  12  except  for  those  events  located  deeper  than  100  km  by  USCGS.  Several  of  the 
earthquakes  on  the  explosion  side  of  the  decision  line  on  Fig.  12  were  deep  and  do  not 
appear  on  Fig.  15.  This  is  quite  reasonable  since  P  waves  from  deep  events  travel 
through  less  low  Q  material  to  reach  the  receiving  station.  It  may  also  be  possible 
that  the  deep  events  actually  generate  more  high  frequency  signals  than  do  normal 
depth  events.  In  any  case  it  is  clear  that  spectral  ratio  does  not  operate  well  for  deep 
events.  Somehow  deep  events  must  be  identified  as  such  and  excluded  from  further 

4 

consideration.  This  is  consistent  with  earlier  results. 

In  Section  VI  it  was  suggested  that  M  —  m,  discrimination  might  be  improved 
by  the  use  of  USCGS  m^  rather  than  LASA  m^.  Our  study  of  the  data  has  not  unearthed 
any  such  preference  in  the  case  of  spectral  ratio.  Either  the  LASA  or  the  USCGS 
seem  to  be  equally  effective  for  discrimination. 

In  addition  to  the  deep  events,  there  are  still  others  which  should  be  excluded 
from  spectral  ratio  tests.  These  are  events  for  which  the  spectral  ratio  is  too  extremely 
distorted  by  the  background  noise  to  be  meaningful.  This  problem  is  considered  in 
Section  VIII.  It  will  be  shown  in  particular  that  the  two  presumed  explosions  incor¬ 
rectly  classified  on  Fig.  15  should  be  discounted  since  the  SNR  is  too  low  for  those 
events. 
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VIII.  MODIFIED  SPECTRAL  RATIO  DISCRIMINATION 


The  spectral  ratio  measured  for  any  event  is  in  fact  a  measure  for  the  sum  of 
signal  and  background  noise.  Ideally  there  would  be  no  background  noise.  This  con¬ 
dition  is  essentially  obtained  for  very  large  events.  However,  it  is  clear  that  the 
measured  spectral  ratio  must  be  increasingly  affected  by  noise  as  events  get  smaller. 
In  the  limit  the  measured  spectral  ratio  is  a  measure  of  noise  properties  and  has 

nothing  at  all  to  do  with  the  signal.  This  signal-to-noise  ratio  problem  has  been  pre- 

4 

viously  noted  by  Kelly.  He  estimated  that  spectral  ratio  could  be  reliably  applied  to 

events  with  m,  §  4.  5.  The  effects  of  noise  contamination  have  now  been  studied  in 
b 

greater  detail  and  results  described  below.  One  fact  which  has  become  apparent  is 
that  the  validity  of  individual  measurements  can  be  determined  by  considering  the 
noise  preceding  the  event. 

The  short -period  beams  for  all  events  are  recorded  on  the  composite  event 
tape.  In  most  cases  about  30  seconds  of  data  is  included  preceding  the  event  arrival. 

It  has  thus  been  possible  to  obtain  two  10-second  noise  intervals  immediately  pre¬ 
ceding  each  event  as  well  as  a  10-second  interval  containing  the  event.  The  intervals 
do  not  overlap.  These  three  intervals  have  been  used  to  study  the  effect  of  noise 
upon  spectral  ratio. 

Recall  that  and  were  previously  taken  to  be  the  high-and  low-frequency 
band  measurements  whose  ratio  is  defined  to  be  spectral  ratio.  Now  let  V^, 
and  Vpjg  be  the  value  of  the  10-second  measured  on  a  beam  for  the  first  noise  inter 
val,  the  second  noise  interval,  and  the  signal  interval,  respectively.  Let  V^> 
and  be  similarly  defined  but  in  terms  of  V^.  Thus,  for  example,  V^/V^  and 
Vh2^L2  wou^  be  spectral  ratio  measurements  made  on  noise  only. 
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We  can  now  define  the  signal-to-noise  factors  =  V^g/V^  and  p^  =  ^H3^H2 
which  are  indicative  of  the  extent  to  which  noise  might  be  effecting  V^g.  The  similarly 
defined  quantities  pj^  and  p^g  can  be  used  for  the  low-frequency  band.  Figures  16 
and  17  show  plots  of  these  quantities  as  a  function  of  LASA  magnitude.  Each  event  is 
represented  by  a  vertical  bar  with  extremes  given  by  (p^,  or  (P^i*  Pl*?)*  ^ie 
relatively  small  variations  in  measured  p  for  a  given  event  (typically  less  than  50%)  is 
accounted  for  by  the  stability  of  the  measurements  on  noise  preceding  the  event. 

Figure  18  shows  histograms  for  some  measurements  of  an<3  V^/V^g  which 

verify  this  stability. 

The  signal-to-noise  factor  is  clearly  not  the  same  as  signal-to-noise  ratio. 

2 

Figure  19  is  a  plot  of  10  log  (p  —  1)  vs  p,  which  crudely  relates  the  noise  factor  p  to 
the  signal-to-noise  ratio  (SNR).  The  relationship  was  obtained  by  assuming  that  both 
the  noise  and  signal  spectra  are  white  noise  random  processes  with  frequency  as  the 
parameter.  The  figure  actually  relates  the  ratio  of  expected  values  of  the  V's  to  the 
signal-to-noise  ratio.  Although  the  relationship  is  admittedly  crude,  it  is  convenient 
to  be  able  to  relate  p  to  a  more  common  parameter  such  as  SNR.  For  example,  p  =  3 
is  about  the  same  as  9  db  SNR  and  p  =  1.  4  is  the  same  as  about  0  db  SNR. 

Since  SNR  can  seriously  affect  spectral  ratio  it  is  necessary  to  make  decisions 
concerning  the  validity  of  spectral  ratio  measurements.  There  are  four  distinct  SNR 
situations.  First,  if  SNR  is  large  in  both  the  higlrand  low-frequency  bands,  then  the 
spectral  ratio  measurement  is  valid  and  truly  reflects  the  event  characteristics. 

Second,  if  the  SNR  is  small  in  both  bands,  then  the  measurement  of  spectral  ratio  can 
lead  to  incorrect  conclusions  since  the  noise  properties  determine  the  spectral  ratio. 
Third,  if  the  SNR  is  poor  only  in  the  high-frequency  band,  then  the  measured  spectral 
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ratio  will  be  increased  from  the  true  value  by  the  noise.  Thus  if  the  event  appears  in 
the  earthquake  region  with  respect  to  spectral  ratio,  it  would  have  been  in  that  region 
even  in  the  absence  of  noise.  However,  if  the  event  appears  to  be  in  the  explosion 
region,  it  means  no  decision  should  be  made.  The  fourth  situation  is  that  in  which  the 
SNR  is  poor  only  on  the  low-frequency  band.  In  this  case  the  noise  tends  to  reduce  the 
measured  spectral  ratio.  Events  which  appear  in  the  explosion  region  should  be 
accepted  since  the  noise  could  only  depress  the  spectral  ratio.  No  decision  can  be 
made  for  events  which  appear  to  be  in  the  earthquake  region. 

We  have  considered  SNR  large  enough  if  it  is  greater  than  9  db  (p  §  3.0). 

Table  I  summarizes  the  decision  rules  which  take  SNR  into  consideration  if  p  =  3. 0  is 
used  as  a  threshold  value  on  the  signal-to-noise  factor. 

Figure  20  shows  spectral  ratio  vs  LASA  for  the  events  on  Fig.  15  which 
remain  after  the  SNR  test  is  applied.  That  is,  events  which  appear  on  Fig.  15  but  not 
on  Fig.  20  have  SNR  values  which  preclude  making  a  decision  based  upon  spectral 
ratio.  Events  on  Fig.  20  indicate  the  true  power  of  spectral  ratio  as  a  discriminant 
since  the  spectral  ratio  for  those  events  is  only  modestly  effected  by  noise.  Events 
indicated  by  numbers  and  letters  on  Fig.  20  are  discussed  further  in  Section  XII. 

Since  spectral  ratio  cannot  be  applied  to  all  events  independent  of  SNR  it  is 
important  to  know,  as  a  function  of  m^,  what  is  the  probability  that  it  can  be  applied. 
Figure  21  shows  this  probability  under  the  condition  that  the  event  is  an  earthquake. 
The  data  points  on  the  figure  are  the  fractions  of  events  in  0. 1  magnitude  intervals  on 
Figure  15  which  also  appear  on  Fig.  20.  The  probability  can  be  interpreted  as  the 
probability  that  a  good  indication  of  event  type  will  be  generated  by  spectral  ratio  if 
noise  free  spectral  ratio  is  a  good  discriminant. 
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TABLE  1 


Modified  Spectral  Ratio  Decision  Rule 


0H>3’  t>3 

Accept  spectral  ratio  measurement  p, 
and  discriminate  according  to  the  de¬ 
cision  line  shown  on  Fig.  12. 

Ph  >  3’  Pl  <  ^ 

If  p  indicates  explosion,  then  decide 
explosion.  If  p  indicates  earthquake, 
make  no  decision. 

Ph<3,  pL>3 

If  p  indicates  earthquake,  then  decide 
earthquake.  If  p  indicates  explosion, 
make  no  decision. 

PH<  3’  PL<  3 

Make  no  decision. 
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It  was  stated  in  Section  VII  that  the  10-second  spectral  ratio  was  at  least  as 
successful  as  the  20-second  spectral  ratio  as  a  discriminant.  Figure  22  can  be  com¬ 
pared  with  Fig.  20  to  verify  this  assertion.  The  plot  of  20-second  spectral  ratio  vs 
contains  all  the  same  points  as  Fig.  20.  The  comparison  is  made  on  this  data  set  since 
spurious  noise  effects  have  been  minimized  as  much  as  possible. 
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IX. 


LASA  PERIOD  AS  A  DISCRIMINANT 


It  seems  that,  in  the  absence  of  noise,  the  short -period  spectra  of  events  can 
be  used  as  a  discriminant.  However,  it  was  demonstrated  in  the  previous  section  that 
noise  can  seriously  effect  the  use  of  spectral  information  at  small  magnitudes.  This 
motivates  die  consideration  of  a  spectrum  related  discriminant  which  might  be  less 
effected  by  noise.  Dominant  period  is  a  candidate  for  such  a  discriminant. 

The  dominant  period  of  an  event  is  generally  defined  in  terms  of  the  seismogram 
itself  although  one  normally  thinks  of  it  as  a  spectral  property.  The  relationship  be¬ 
tween  spectrum  and  dominant  period  is  at  best  crude  and  will  not  be  discussed  here. 

We  have  taken  dominant  period  to  be  twice  the  time  between  the  positive  and  negative 
parts  of  the  largest  peak-to-peak  excursions  of  LASA  beams  which  have  been  bandpass 
filtered  from  0.  6  to  2.  0  Hz.  SPLP  generates  the  dominant  period  and  is  not  subject  to 
human  subjectivity.  Due  to  the  measurement  procedure  and  the  LASA  sampling  rate 
all  dominant  periods  were  reported  as  multiples  of  0. 1  seconds. 

Figure  23  shows  dominant  period  v£  LASA  magnitude  for  all  events  in  our 
population  except  those  reported  at  depths  greater  than  100  km  by  USCGS.  One  pos¬ 
sible  decision  line  has  been  drawn  which  contains  all  explosions  in  the  lower  right  sec- 

4 

tion  of  the  plot.  This  data  is  similar  to  that  presented  previously  by  Kelly  but  is 
more  complete  and  has  had  some  of  the  human  factors  removed. 

If  Fig.  23  is  accepted  at  face  value,  it  suggests  that  all  explosions,  even  down 
to  below  magnitude  4.  0,  can  be  correctly  indicated  but  that  only  53%  of  the  earthquakes 
will  be  correctly  identified.  However,  it  also  suggests  that  the  quality  of  a  decision 
based  upon  period  might  strongly  depend  upon  magnitude.  If  attention  is  directed  to 
the  important  range  below  magnitude  4.  5,  we  observe  that  all  explosions  are  correctly 
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identified  and  72%  of  the  earthquakes  are  correctly  identified.  If  such  performance  is 
actually  possible,  then  dominant  period  can  be  a  valuable  discriminant  in  the  low  mag¬ 
nitude  range. 

Before  terminating  this  section,  it  is  important  to  note  that  dominant  period  as 
well  as  spectral  ratio  may  be  highly  dependent  upon  the  receiver  location  as  well  as 
the  source  regions.  The  properties  of  spectra  in  the  low-and  high-frequency  bands  as 
well  as  dominant  period  may  be  very  different  for  another  receiver  even  for  the  same 
source  regions. 
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X. 


LASA  COMPLEXITY  AS  A  DISCRIMINANT 


The  value  of  complexity  as  a  discriminant  has  been  quite  extensively  discussed 
4 

by  Kelly.  He  concluded  that  complexity  is  of  limited  usefulness.  Our  data  has  sub¬ 
stantiated  and  strengthened  that  conclusion. 

Figure  24  is  a  plot  of  complexity  vs  magnitude  for  all  events  except  those 
located  at  depths  greater  than  100  km  by  USCGS.  The  data  contains  more  large  earth¬ 
quakes  and  small  explosions  than  have  been  previously  available.  It  is  clear  that  even 
for  large  events  (say  >  4.  8),  for  which  there  are  no  noise  contamination  difficulties, 
complexity  is  not  a  powerful  discriminant.  A  decision  line  is  shown  drawn  at  3.  0  so 
that  all  explosions  are  properly  identified.  Complexity  then  correctly  identifies  only 
53%  of  the  large  earthquakes  and  only  47%  of  all  earthquakes.  If  the  decision  line  is 
brought  down  to  a  complexity  of  2.  0,  the  percent  of  correctly  identified  earthquakes 
increases  to  73%.  However,  two  explosions  are  then  moved  well  into  the  earthquake 
region.  One  can  probably  ignore  the  small  explosion  since  the  background  noise  is 
forcing  the  complexity  to  be  large.  The  large  event  cannot  be  ignored  and,  as  previ¬ 
ously  pointed  out  by  Kelly,  forces  retention  of  the  higher  decision  level  which  results 
in  poorer  performance. 
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XI.  DEPTH  PHASES  AT  LAS  A 


Earthquake  depth  is  an  extremely  important  parameter  in  any  discrimination 
study.  For  example,  if  depth  is  firmly  established  to  be  more  than  just  a  few  kilo¬ 
meters,  an  event  can  obviously  be  identified  as  an  earthquake.  In  addition,  certain 
other  discriminants,  such  as  complexity,  dominant  period,  spectral  ratio  or  the 
M  —  m,  relationship,  are  generally  believed  to  be  a  function  of  depth.  Indeed,  deep 
events  ( >  100  km)  often  generate  very  small  surface  waves  and  often  have  teleseismic 
P  waves  very  much  like  those  of  underground  explosions. 

All  events  in  our  population  have  been  examined  by  an  analyst  for  possible  depth 
phases  (sP  or  pP).  The  analyst  operated  as  follows.  First  the  SPLP  beam  waveforms 
were  examined  for  any  possible  arrivals  after  the  initial  P  wave  onset.  Picks  corre¬ 
sponding  to  PcP  times  were  discarded.  The  time  corresponding  to  the  largest  ampli¬ 
tude  of  any  remaining  picks  was  then  accepted  as  a  candidate  for  sP  or  pP.  The  anal¬ 
yst  accepted  picks  associated  with  amplitudes  less  than  25%  of  the  P  amplitude  only  if 
he  felt  the  arrival  very  clearly  rose  above  the  coda  or  noise  level. 

The  above  procedure  was  scrupulously  executed  for  presumed  explosions  as 
well  as  for  earthquakes.  As  a  result  erroneous  depth  phases  were  attributed  to  three 
explosions  in  the  population.  Two  of  these  were  complex  Novaya  Zemlya  events  with 
many  apparent  secondary  arrivals,  some  as  large  as  0.8  of  the  initial  P  amplitude. 

The  third  event  was  quite  simple  and  is  s  hown  as  event  2  in  Fig.  26. 

The  danger  of  accepting  apparent  depth  phases  at  LASA  as  unequivocal  evidence  of 
depth  for  discrimination  is  obvious.  Corroborating  evidence  from  other  sources  should 
be  required  in  order  to  avoid  the  incorrect  classification  of  an  explosion  as  an  earth¬ 
quake. 
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The  analyst  accepted  depth  phases  for  58%  of  all  the  earthquakes  in  our  data 
base.  This  is  just  slightly  less  than  the  acceptance  rate  for  the  55%  of  all  earthquakes 
in  the  data  base  which  were  also  reported  by  USCGS.  For  that  subset  of  earthquakes 
depth  phases  were  attributed  to  63%  of  the  events  by  the  analyst. 

If  events  restrained  to  33  km  by  USCGS  are  excluded  from  consideration  there 
are  41  earthquakes  in  our  population  to  which  we  attributed  depth  phases  and  which 
were  reported  by  USCGS.  Figure  25  shows  the  degree  of  agreement  between  USCGS 
depths  and  the  depths  implied  by  the  possible  depth  phases  observed  on  the  SPLP  beams. 
Each  event  is  represented  by  (1)  a  vertical  line  indicating  the  USCGS  confidence  inter¬ 
val  and  (2)  a  solid  dot  indicating  the  SPLP  depth  which  is  obtained  by  assuming  the 
phase  is  pP.  In  cases  where  a  depth  computed  by  assuming  the  phase  to  be  sP  gives 
better  agreement  with  USCGS  the  sP  depth  has  been  indicated  by  an  open  dot.  In  cases 
where  there  might  be  confusion  lines  connecting  dots  and  vertical  lines  have  been  in¬ 
cluded  to  indicate  that  they  should  be  associated  with  the  same  event. 

Given  the  available  data  shown  on  Fig.  25  it  is  at  least  difficult  to  precisely 

evaluate  the  success  of  the  depth  determinations  made  using  SPLP.  There  are  only 

two  cases  (LASA  pP  depth  =  29  km,  USCGS  depth  =  87  km;  LASA  pP  depth  =  80;  USCGS 

depth  =  10)  in  which  there  is  virtually  no  doubt  that  the  LASA  phase  was  neither  pP  nor 

sP.  There  are  about  30  cases  in  which  it  is  most  reasonable  to  assume  that  the  SPLP 

phase  was  either  pP  or  sP.  The  remaining  events  are  ambiguous  at  best.  Thus,  for 

earthquakes,  between  70%  and  95%  of  the  possible  depth  phases  picked  on  SPLP  beams 

corresponded  to  either  sP  orpP.  This  implies  that  sP  or  pP  was  correctly  picked  for 
between  42%  and  57%  of  all  the  earthquakes  in  our  population. 
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Observe  that  the  decision  based  upon  a  possible  depth  phase  is  the  only  one 
which  is  in  error  by  incorrectly  identifying  a  presumed  explosion  as  an  earthquake. 

If  such  an  error  is  very  serious,  then  discrimination  based  upon  depth  phases  is  not 
acceptable  without  much  more  verification  than  is  possible  at  a  single  site.  However, 
the  fact  that  all  but  one  of  the  earthquakes  which  were  problems  for  spectral  ratio 
appear  to  have  depth  phases  is  suggestive  of  the  worth  of  depth  phases  if  they  can  be 
adequately  verified.  Because  of  this  problem  of  the  false  classification  of  an  explosion 
we  will  not  consider  possible  depth  phases  as  a  discriminant  in  the  remainder  of  this 
section. 

It  is  obvious  from  Table  III  that  not  all  discriminants  yield  the  same  decision 
for  a  given  event.  Even  leaving  possible  depth  phases  out  of  consideration  there  are 
clearly  disagreements  for  events  C,  E,  G,  I.  For  example,  events  C  and  G  are  iden¬ 
tified  as  earthquakes  by  Mg  —  and  as  explosions  by  spectral  ratio.  Event  E  is 
identified  clearly  as  an  earthquake  by  first  motion  but  is  clearly  indicated  to  bean 
explosion  by  the  LASA  Mg  —  criteria.  Event  I  is  apparently  incorrectly  identified 
by  spectral  ratio  but  is  correctly  identified  by  complexity.  In  addition,  the  first 
motion  may  be  dilatational  although  we  have  not  indicated  it  as  such  in  table  III.  We 
hope  that  these  various  disagreements  adequately  demonstrate  the  multivariat  nature 
of  the  discrimination  problem. 
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TABLE  II 


Symbol  of 
Figs.  20,22 

mb 

M  * 
s 

Complexity 

Period 

First 
Motion  ** 

Time  to 
Possible  Phase 

Spectral 

Ratio 

A 

6.  24 

4.  76 

3.1 

0.8 

- 

8 

0.  39 

B 

5.  75 

4.45 

5.0 

0.8 

+ 

13 

0.  45 

C 

5.81 

4.20 

2.6 

0.7 

+ 

12 

0.  79 

D 

w 

5.87 

5.63 

5.  7 

1.  0 

+ 

13 

0.39 

0) 

17  M 

E  g 

5.70 

3.21 

1.6 

0.9 

— 

10 

0.  40 

a* 

4.43 

<3.  53 

3.3 

0.8 

+ 

11 

0.  76 

G 

4.  82 

3.09 

2.0 

0.8 

+ 

5 

0.86 

H 

5.44 

3.64 

1.4 

0.9 

- 

0.  45 

I 

4.  65 

<4.94 

3.3 

0.8 

+ 

23 

1.00 

J 

4.  91 

7.1 

1.0 

+ 

11 

0.  63 

'll 

4.  69 

<2.  75 

1.59 

0.8 

+ 

0.  77 

2s§ 
w  2, 

g  % 

5.87 

3.57 

1.38 

0.9 

+ 

20 

0.  42 

3 

4.  95 

<3.  30 

0.  89 

0.8 

+ 

— 

0.89 

*  <  indicates  an  upper  limit  on  M  imposed  by  noise  or  an  interfering  event 

s 

**  First  motion  is  indicated  +  if  it  is  +  or  highly  ambiguous. 


XII.  DETAILED  DATA  FOR  SOME  SPECIFIC  EVENTS 


It  is  the  behavior  of  individual  discriminants  which  has  been  considered  in  the 
preceding  sections.  In  this  section  the  data  from  some  specific  events  are  used  to 
give  some  indication  that  discrimination  is  not  just  a  one -dimensional  problem.  The 
interaction  and  statistical  dependence  between  various  discriminants  is  of  great  impor¬ 
tance.  Although  we  have  not  investigated  this  area  in  depth  it  is  hoped  that  our  con¬ 
sideration  of  several  discriminants  for  a  few  events  will  serve  to  indicate  how  impor¬ 
tant  it  is. 

Three  presumed  explosions  and  10  earthquakes  were  singled  out  on  Fig.  20 
because  they  were  either  incorrectly  identified  by  spectral  ratio  or  were  the  most 
marginal  events.  Figure  26  gives  waveform  information  for  those  events.  The 
SPLP  data  is  summarized  in  Table  II.  Since  first  motion  is  obviously  important  for 
these  events  the  table  includes  a  column  for  the  polarity  of  first  motion  as  well  as  for 
the  discriminants  discussed  in  earlier  sections. 

Table  III  shows  decisions  which  might  be  made  using  the  data  in  Table  II.  A 
decision  in  favor  of  the  earthquake  hypothesis  is  indicated  by  "quake.  "  A  decision  in 
favor  of  the  explosion  hypothesis  is  indicated  by  "exp.  "  A  blank  appears  when,  in  our 
opinion,  the  value  of  the  discriminant  does  not  suggest  one  hypothesis  much  more 
strongly  than  it  suggests  the  other. 

The  period  column  of  Table  III  is  completely  blank.  If  the  period  of  any  of  the 
events  had  been  sufficiently  long  (see  Fig.  23)  then  an  earthquake  decision  would  have 
been  made.  However,  a  short -period,  except  at  magnitudes  lower  than  that  of  any  of 
the  13  events  under  consideration,  does  not  significantly  reduce  the  probability  that 
the  event  is  an  earthquake.  Thus  we  have  decided  to  indicate  no  decisions  based  upon 
period. 
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TABLE  III 


Event 

Symbol 

M  — m, 
s  b 

Complexity 

Period 

Spectral 

Ratio 

First 

Motion 

Possible 

Depth 

Phase 

A 

Quake 

Quake 

Quake 

B 

Quake 

Quake 

Quake 

C 

Quake 

Exp 

Quake 

D 

Quake 

Quake 

Quake 

E 

Exp 

Quake 

Quake 

F 

Quake 

Quake 

G 

Quake 

Exp 

Quake 

H 

Quake 

Quake 

I 

Quake 

Exp 

Quake 

J 

Quake 

Quake 

1 

2 

Exp 

Quake 

3 

Exp 

XIII.  SUMMARY 


The  primary  LASA  digital  recordings  for  nearly  200  events  have  been  analyzed 
to  obtain  refined  data  for  a  variety  of  discrimination  experiments.  Some  of  the  dis¬ 
crimination  studies  have  been  completed  and  are  reported  in  detail  in  the  body  of  this 
report.  The  report  has  been  primarily  limited  to  the  evaluation  of  individual  dis¬ 
criminants  measured  at  a  single  site.  In  the  most  general  terms  our  data  corroborates 
the  discrimination  performance  of  M  —  m,  ,  spectral  ratio,  period,  and  complexity 

S  D 

4  5 

which  has  been  previously  reported  by  Kelly  and  Capon,  et  al. 

The  ability  of  LASA  to  discriminate  on  the  basis  of  Mg  —  has  been  evaluated 
for  those  events  which  did  not  have  other  events  interfering  with  the  LP  data  and  which 
were  less  than  about  100  km  in  depth.  One  Kuriles  earthquake  appeared  as  an  explosion 
unless  the  USCGS  value  of  was  used.  Ignoring  this  one  anomalous  event  LASA 
discriminated  perfectly  for  events  with  m^  1  4.  8  in  all  cases  where  surface  waves  were 
detected.  Exclusive  of  events  obscured  by  other  events  or  for  which  no  usable  tape 
existed,  it  was  possible  to  observe  surface  waves  for  all  earthquakes  with  m^  5  4.  8. 

All  presumed  explosions  with  m^  5  5.1  could  also  be  identified  using  their  My  value 
or  the  bound  imposed  by  the  background  noise.  Below  m^  5  5.1  no  events  were  un¬ 
equivocally  identified  as  explosions.  However,  the  probability  that  an  earthquake  could 
be  identified  as  such  decreased  to  zero  only  gradually  as  was  reduced  from  4. 8  to 
about  4. 0. 

The  spectral  ratio  criteria  has  been  modified  to  utilize  estimates  of  the  signal - 
to-noise  ratio  in  the  bands  used  to  construct  the  spectral  ratio.  The  modification  intro¬ 
duces  the  option  to  make  no  decision.  Hie  probability  that  no  decision  is  made  is  zero 
for  magnitudes  above  about  4.  8  and  unity  below  about  4. 0.  Once  a  decision  is  made. 
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the  probability  that  it  is  correct,  given  that  events  deeper  than  100  km  are  excluded, 
is  very  high. 

Period,  plotted  against  m^,  may  be  a  discriminant  which  will  be  valuable  at  low 
magnitudes  (say  below  4.  5).  It  has  not  been  possible  to  unequivocally  identify  events 
as  explosions  using  period.  However,  many  earthquakes  with  magnitudes  below  4.  5 
can  be  identified  correctly  using  the  period  data. 

Event  complexity  measured  at  LASA  appears  to  be  a  relatively  weak  discriminant. 
Noise  effects  negate  its  value  for  small  events.  Large  events  tend  to  be  identified  by 
Mg  —  m^,  spectral  ratio,  and  depth  phases.  Qf  course  all  non-trivial  discriminants 
should  play  a  role  in  multi-variate  or  multi-station  discrimination.  Complexity  may 
occasionally  be  of  more  value  in  that  context. 

Possible  depth  phase  picks  have  been  made  for  about  60%  of  the  earthquakes  in 
our  population.  A  comparison  with  USCGS  depths  has  indicated  that  between  70%  and 
95%  of  the  picks  are  valid  depth  phases.  Unfortunately,  depth  phases  were  also  picked 
for  three  presumed  explosions. 

All  of  our  data  is  from  a  single  array  station  monitoring  a  limited  part  of  the 
world.  The  evaluation  of  discriminants  is  strictly  valid  only  for  those  conditions.  Of 
course  there  are  general  physical  explanations  why  the  various  discriminants  should 
work  to  some  extent.  However,  the  success  may  vary  with  site  and  it  may  be  neces¬ 
sary  to  significantly  modify  details  to  achieve  successful  operation  at  different  sites. 

For  example,  it  is  anticipated  that  differences  in  earthquake  and  explosion  short- 
period  spectra  will  appear  at  other  sites  such  as  NORSAR,  but  that  spectral  ratio  as 
presently  defined  may  not  make  the  most  effective  use  of  those  differences. 
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APPENDIX  A 


FREQUENCY  DOMAIN  FILTERING  OPERATIONS 


All  bandpass  and  chirp  filtering  operations  mentioned  in  this  report  were 
accomplished  using  Discrete  Fourier  Transforms  (DFT).  In  this  Section  we  briefly 
describe  the  DFT  and  the  filtering  procedure  which  has  been  used.  Extensive  exposi¬ 
tions  of  the  DFT  and  its  properties  are  readily  available  in  the  literature. 

Let  s(n),  n  =  0,  1, ...  ,N  —  1  be  samples  of  a  seismogram  which  we  wish  to 
filter.  The  DFT  of  this  sequence  is  defined  as 


N-l 


s(n)  e'i2n(nk/N> 


S(k)  =  E 
n=0 


It  is  assumed  that  k  takes  on  integer  values  only.  The  function  S(k)  is  periodic  with 
period  N.  That  is, 


S(—  k)  =  S(N-k)  . 


If  the  sequence  s(n)  is  periodically  extended  outside  the  interval  [0,  N  —  1]  by  the 
relation 

s  ( —  n)  =  s(N  —  n) 

then  s(n)  is  obtained  from  S(k)  by  inverse  transformation  according  to 


1  N'1 

s(n)  =  i  Z  S(k)  e 
iN  k=0 


i2rr(nk/N) 
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Suppose  that  the  time  between  s(n)  and  s(n  +  1)  is  6  seconds.  That  is,  s(n), 
n  =  0, . . .  ,N  —  1  can  be  considered  as  a  sequence  of  samples  taken  every  6  seconds 
from  a  seismogram.  The  Nyquist  frequency  for  this  case  is  1/26  cycles  per  second. 
Using  this,  one  can  associate  the  integer  argument  of  a  DFT  on  the  interval  [0,N  —  1] 
with  real  time  positive  and  negative  frequencies.  Specifically 


for  k  =  0,  l,...,N/2  — 1  and 


-  N-k 
-fk  ”  ~"N§ 

for  k  =  N/2  +  1,. . .  ,N  —  1. 

A  simple  but  effective  bandpass  filter  can  be  applied  to  s(n)  by  setting  all  S(k) 
to  zero  for  k  which  correspond  to  frequencies  outside  of  the  band.  For  example,  let 
fpj  and  f-^  be  the  high  and  low  edge  of  the  band.  Define  k^  =  N6f^  and  k^  =  N6f^.  We 
assume  that  k^  and  k^  are  positive,  are  founded  to  the  nearest  integers,  and  that 
<<:  N/2*  The  bandpassed  signal  is  now  obtained  as 

b(n)  =  i  Z  S(k)  el2TT(nk/N>  +  S(N  -  k)  el2Tl[n(N  ~  k)/N]  . 

R=kL 

Because  s(n)  is  real  it  follows  that  b(n)  is  real  also. 

Except  for  modifications  due  to  the  periodicity  of  S  and  periodic  extension  of  s 
the  DFT  behaves  in  many  ways  like  the  ordinary  Fourier  transform.  In  particular, 
multiplication  of  functions  in  the  frequency  domain  corresponds  to  convolution  in  the 
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time  domain.  We  shall  refer  to  the  inverse  transform  of  the  frequency  characteristic 
of  any  filter  as  the  impulse  response.  We  have  only  considered  filters  with  frequency 
functions  H(k)  which  satisfy  H(k)  =  H(N-k)  so  their  impulse  responses  are  real.  Thus 
b(n)  above  can  be  obtained  by  convolving  the  periodically  extended  s(n)  with  the  impulse 
response 


h(n)  =  ‘  EH  ei2n<nk/N>  +  ei2T,["<N "  W™ 
k=KL 


The  function  h(n)  is  shown  in  Fig.  27  for  N  -  2048,  6  =  1.0,  f ^  =  0.  055  and  f^  =  0.  025. 

These  parameters  were  used  when  filtering  long-period  data.  Because  h(n)  introduces 

no  phase  distortion,  it  is  symmetrical  about  n  =  0  as  well  as  periodic  with  period  N. 

The  impulse  response  is  essentially  zero  outside  of  the  time  interval  shown. 

13 

Chirp  filters  are  most  often  considered  in  the  time  domain  in  terms  of  their 
impulse  response.  However,  the  fast  algorithms  for  DFT  make  it  desirable  to  perform 
chirp  filtering  operations  in  the  frequency  domain.  We  have  done  chirp  filtering  using 
frequency  response  functions 


G(k)  = 


"  ei2TT(C/N)  <k  -  ko)^  u  ^  g  k  s  kw 
,  e-i2TT(C/N)(N-k-k ^HN-vstSN-t, 

o  Lf 

0  ,  otherwise 


Such  a  frequency  function  will  pass  energy  in  the  band  implied  by  k^  and  k^  while 
rejecting  energy  outside  the  band.  The  phase  shift  through  the  band  is  quadratic  with 
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frequency.  There  will  be  zero  phase  shift,  resulting  in  zero  delay,  of  the  frequency 
component  fQ  =  N6kQ.  Such  a  frequency  function  does  not  exactly  correspond  to  the 
usual  chirp  filter  impulse  response,  which  is  a  burst  of  sine  wave  with  linearly  changing 
frequency.  This  is  not  of  any  real  importance  as  long  as  C  can  be  adjusted  to  achieve 
from  G  a  filter  reasonably  matched  to  the  waveforms  of  interest  (Dispersed  Rayleigh 
wave  trains). 

Figure  28  shows  impulse  responses  corresponding  to  two  typical  values  of  C 
which  have  been  used.  The  zero  delay  frequency  used  is  0.  05  Hz  and  N  =  2048  with 
6=1.0  second.  The  upper  and  lower  bandpass  frequencies  are  0.  055  Hz  and  0.  025  Hz 
respectively.  Outside  the  time  interval  shown  the  waveforms  have  even  smaller  ampli¬ 
tudes  than  at  the  edges.  The  amplitude  outside  the  interval  shown  can  be  considered 
to  be  zero. 
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APPENDIX  B 


Each  record  of  the  composite  event  tape  contains  exactly  4131  four-byte  words. 
The  records  can  be  directly  read  by  a  binary  read  in  a  Fortran  IV  program  operating 
in  an  IBM  360  system.  Each  word  of  the  record  is  properly  interpreted  in  the  program 
as  either  an  integer  or  a  floating  point  number  as  indicated  in  the  following  tabular 
description  of  the  records  on  the  tape.  All  data  is  as  determined  by  LASA  and  SPLP. 


WORD  TYPE  (INTEGER  OR 

FLOATING  POINT) 

DESCRIPTION 

1  I 

Arbitrarily  assigned  event  number.  Numbers 
in  the  1500  to  1999  range  are  possible  explo¬ 
sions 

2-7  I 

Hour,  minute,  second,  tenths  of  seconds, 
day  number,  and  year  for  event  arrival  at 
the  center  of  LASA.  For  impulsive  events 
this  is  the  time  of  peak  excursion.  For 
others,  it  is  an  approximate  onset  time. 

8  F 

Zero-to-peak  amplitude  of  P  wave  (m|_i) 

9  F 

Period  of  P  wave  (seconds) 

10  F 

Magnitude 

11  F 

P-wave  azimuth 

12  F 

P-wave  horizontal  phase  velocity 
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WORD 


TYPE  (INTEGER  OR 
FLOATING  POINT) 


DESCRIPTION 


13  F 

14  F 


15  F 

16  I 


17  F 

18  I 


19 


F 


20 


I 


21  F 

22  I 


23 


F 


Latitude  (degrees  N) 

Longitude  (degrees  E) 

Distance  (degrees) 

This  word  is  set  to  unity  if  the  peak  excursion 
of  the  SP  P-wave  is  more  than  2.  5  seconds 
after  onset;  otherwise  this  word  is  zero. 

This  is  the  surface  factor 

This  word  is  zero  if  a  good  surface  wave 
wave  detection  was  achieved.  It  is  unity  if 
normal  LP  background  noise  obscured  the 
event.  It  is  two  if  surface  waves  from  an¬ 
other  event  obscured  the  event  of  interest. 

Any  other  value  implies  that  either  no  LP 
data  was  available  or  equipment  errors 
made  it  unusable. 

This  word  is  zero  unless  a  possible  depth 
phase  was  picked.  In  that  case  it  equals  the 
moveout  to  the  possible  phase  (seconds). 

Number  of  0.  05  second  digital  samples  before 
the  time  given  in  words  2  —  7  to  the  start  of 
the  complexity  measurement  interval. 

Complexity 

Number  of  samples  before  the  time  given  in 
words  2  —  7  to  the  beginning  of  the  interval 
used  for  the  measurement  of  spectra  and 
spectral  ratio. 

Spectral  ratio  measured  using  a  10-second 
interval  of  data. 


45 


WORD 


TYPE  (INTEGER  OR 
FLOATING  POINT) 


DESCRIPTION 


24 


F 


Spectral  ratio  measured  using  a  20-second 
interval  of  data. 


25-30  I 

31  I 

32-2079  F 

2080  I 

2081-3105  F 

3106  I 

3107-4131  F 


Hour,  minute,  second,  tenth  of  seconds,  day 
number,  and  year  for  the  first  sample  of  the 
short -period  beam. 

Number  of  0.  05  second  samples  in  the  beam 
(NBM) 

The  best  unfiltered  short -period  beam.  Data 
after  the  first  NBM  samples  are  set  to  zero. 

One  less  than  the  number  of  0. 05  second 
samples  in  the  interval  transformed  to  obtain 
a  spectrum  for  10  seconds  of  data  (NS PR  10). 

Voltage  spectrum  using  10-second  window  on 
the  data.  Data  after  the  first  NS  PR  10  is  set 
to  zero.  Frequencies  range  from  0.  0  to  10.  0 
Hz  in  increments  of  10.  0/NSPR10. 


Same  as  2080-3105  but  for  a  20-second 
interval  of  data. 
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Fig.  1.  Body  wave  magnitude  distribution  of  earthquakes  in  experiment. 
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Fig.  2.  USCGS  vs  LASA  for  events  with  USCGS  depth  less  than  100  km. 
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ir 
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COMPONENT  OF  LP  DATA  FROM 
UP  TO  21  SITES 


L 


REPEAT  FOR 
EACH  OF  3 
LONG-PERIOD  / 
COMPONENTS 
(2,  N-S,  E-W) 


CHECK  DATA  QUALITY  AND 
REJECT  BAD  DATA 


BANDPASS  FILTER  DATA 
0.025  -0.055  Hz 


FORM  BEAM  STEERED 
TO  NOMINAL  RAYLEIGH 
VELOCITY  OF  3.  7  km/iec 


COMPUTE  HORIZONTAL 
SEISMOGRAMS  RADIAL 
AND  TRANSVERSE  TO 
DIRECTION  OF  ARRIVAL 
OF  EVENT 


APPLY  FAMILY  OF  CHIRP 
FILTERS  TO  LP2  BEAM 
AND  SAVE  THAT  WITH 
LARGEST  ZERO-TO-PEAK 
VALUE 


V 


COMPUTE  SURFACE  WAVE  MAGNITUDE 
FROM  CHIRP  FILTERED  LPZ  BEAM 


Fig.  3.  Sequence  of  operations  in  SPLP  program. 
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LASA 


Fig.  4.  LASA  surface  wave  magnitude  (Ms)  vs  LASA  body 
wave  magnitude  (m^)  for  earthquakes  with  detected  surface 
waves  (events  with  USCGS  depth  k  100  km  excluded). 
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Fig.  5.  Surface  wave  detection  data  for  Sino-Soviet  earthquakes 
recorded  at  LAS  A. 


52 


DETECTION  PROBABILITY  DETECTION  PROBABILITY 


Fig.  6.  Surface  wave  detection 
probability  for  Sino-Soviet  earth¬ 
quakes  (events  with  USCGS  depth 
^100  km  excluded). 


Fig.  7.  Apparent  surface  wave 
magnitude  of  long-period  noise 
(Sino-Soviet  locations  assumed). 


Fig.  8.  Surface  wave  detection 
probability  inferred  from  noise 
observations. 
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Fig.  9.  LASA  Ms  vs  USCGS  mb  for  all  events  with  detected  surface  waves. 
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Fig.  10.  LASA  Ms  vs  LASA  for  events  with  detected  surface  waves 
(events  with  USCGS  depth  £  100  km  excluded). 
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Fig.  11.  LASA  Ms  vs  USCGS  mb  for  detected  events  with  USCGS  depth  <  100  km. 
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Fig.  12.  Spectral  ratio  vs  LASA  for  all  events. 
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Fig.  13.  High-frequency  band  signal  factor  vs  LASA  mb. 
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Fig.  14.  Low-frequency  band  signal  factor  vs  LASA 
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Fig.  15.  Spectral  ratio  vs  LASA  mj,  (events  with  USCGS  depth  ^  100  km  excluded). 
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Fig.  16.  High-frequency  band  signal-to-noise  factor  vs  LASA 
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Fig.  17.  Low-frequency  band  signal-to-noise  factor  vs  LASA  m^. 


HIGH- BAND  NOISE  STABILITY  FACTOR  (VH,/VH2) 


CO 

2 


3 

2 


LOW -BAND  NOISE  STABILITY  FACTOR  (VLJ/VL2) 


Fig.  18.  Stability  study  for  signal-to-noise  factor  test. 


Fig.  19.  Signal-to-noise  ratio  as  a  function  of  signal-to-noise  factor. 
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Fig.  20.  Spectral  ratio  vs  LASA  for  events  admitted  by  signal-to-noise 
ratio  test  (events  with  USCGS  depth  i;  100  km  excluded). 
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Fig.  21.  Decision  probability  for  modified  spectral  ratio  test. 


Fig.  22.  Spectral  ratio  from  20  seconds  of  data  vs  LASA  m^. 
(Same  events  as  Fig.  21.) 
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Fig.  23.  Period  of  0.  6  -  2.  0  Hz  filtered  beam  vs  LASA 
(events  with  USCGS  depth  ^  100  km  excluded). 
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Fig.  24.  Complexity  measured  on  0.  6  -  2. 0  Hz  filtered  beam  vs 
LASA  mj,  (events  with  USCGS  depth  §  100  km  excluded). 
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(a)  Shallow  events. 


Fig.  25.  USCGS  depth  vs  LASA  depth  for  earthquakes  with  possible 
depth  phases  at  LASA. 
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(b)  Deep  events. 


Fig.  25.  Continued. 
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Fig.  26.  Waveforms  for  the  most  ambiguous  events  according  to  the  modified  spectral  ratio  criterion. 
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Fig.  27.  Long  period  band-pass  filter  impulse  response. 


Fig.  28.  Long  period  chirp  filter  impulse  response. 


71 


UNCLASSIFIED 
Security  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 

(Security  ctaaaili  cation  ot  title,  body  ot  ebatract  mnd  Indexing  annotation  muat  be  entered  when  the  overall  report  ia  claaeified) 

1.  ORIGINATING  ACTIVITY  (Corporate  author)  2a.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

Lincoln  Laboratory,  M.I.T.  2b.  group 

None 

3.  REPORT  TITLE 

A  Large-Population  LASA  Discrimination  Experiment 


4.  DESCRIPTIVE  NOTES  (Type  ot  report  and  incluaive  datee) 

Technical  Note 

3.  AUTHOR(S)  (La at  name,  (iret  name,  initial) 


Lacoss,  Richard  T. 


6.  REPORT  DATE 

8  April  1969 


7a.  TOTAL  NO.  OF  PAGES 

76 


7b.  NO.  OF  REFS 

13 


8a.  CONTRACT  OR  GRANT  NO. 

AF  19  (628)-5167 

b.  PROJECT  NO. 

ARP  A  Order  512 


9a.  ORIGINATOR'S  REPORT  NUMBER(S) 

Technical  Note  1969-24 


96.  OTHER  REPORT  NO(S)  (Any  other  ntmbere  that  may  be 
aeaifned  thia  report) 

ESD-TR-69-69 


10.  AVAILABILITY/LIMITATION  NOTICES 


This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 


12.  SPONSORING  MILITARY  ACTIVITY 

Advanced  Research  Projects  Agency , 
Department  of  Defense 

13.  ABSTRACT 


11.  SUPPLEMENTARY  NOTES 

None 


A  computer  program  has  been  written  and  has  been  applied  to  LASA  time  series  data  from  nearly 
200  events  in  order  to  obtain  data  for  discrimination  studies.  In  general,  the  previously  estimated 
performance  of  Ms  —  m,  spectral  ratio,  and  complexity  have  been  corroborated.  Some  of  the  specific 
results  are  the  following.  Except  for  one  anomalous  earthquake,  all  shallow  earthquakes  with  m^  >4.8 
and  explosions  with  m^  >5.1  have  been  correctly  identified  using  Ms  —  m\y  No  events  with  mb  <5.1 
could  be  unequivocally  identified  as  explosions.  However,  the  probability  that  an  earthquake  would  be 
identified  as  such  decreased  to  zero  only  gradually  as  mb  was  reduced  to  about  4.0.  A  modification  of 
spectral  ratio  has  been  made  which  introduced  the  option  to  make  no  decision  concerning  the  nature  of 
an  event  when  the  signal-to-noise  ratio  is  not  sufficiently  large.  The  probability  that  no  decision  will 
be  made  is  zero  for  mb  >4.8  and  unity  for  mb  <4.0.  The  probability  of  correct  identification  is  high 
for  events  which  pass  the  signal-to-noise  ratio  tests.  The  period  of  short  period  P-waves  may  be  a 
valuable  discriminant  at  low  magnitudes.  Many  earthquakes  with  magnitudes  below  4.5  can  be  identified 
as  such  using  the  period  data.  Depth  phase  picks  have  been  made  for  about  60  percent  of  the  earth¬ 
quakes  in  our  population.  About  70  to  95  percent  of  these  picks  correspond  to  valid  depth  phases. 
Unfortunately,  depth  phases  were  also  picked  for  three  presumed  explosions. 


14.  KEY  WORDS 


LASA  earthquakes 

seismic  array  earth  sensors 

seismic  discrimination 


72 


UNCLASSIFIED 


Security  Classification 


4 


